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CN . Some properties of the inverse error function 

^ ■ DieETO Dominici 

I— J 

Abstract. The inverse of the error function, inverf(x), has applications in 
diffusion problems, chemical potentials, ultrasound imaging, etc. We analyze 
' the derivatives inverf (z) , as n — > oo using nested derivatives and 

, a discrete ray method. We obtain a very good approximation of inverf (x) 

• ' through a high-order Taylor expansion around x = 0. Wo give numerical 

' results showing the accuracy of our formulas. 

1. Introduction 

^ ' The error function erf(z), defined by 

o ■ 
m 
(N 

: ' 

■ occurs widely in almost every branch of applied mathematics and mathematical 
physics, e.g., probability and statistics [Wal50], data analysis [Her88], heat con- 

■ duction [Jae46], etc. It plays a fundamental role in asymptotic expansions [01v97] 
I and exponential asymptotics [Ber89]. 

■ Its inverse, which we will denote by inverf (z) , 



z 

erf (2) — —j= \ exp (— i^) dt, 



inverf (z) = erf "'"(z). 



appears in multiple areas of mathematics and the natural sciences. A few examples 
k> \ include concentration-dependent diffusion problems [Phi55], [Sha73], solutions to 

• Einstein's scalar-field equations [LW95], chemical potentials [TM96], the distri- 

5^ I bution of lifetimes in coherent-noise models [WM99] , diffusion rates in tree-ring 

chemistry [BKSH99] and 3D freehand uhrasound imaging [SJEMFAL+03]. 

Although some authors have studied the function inverf (z) (see [DomOSb] 
and references therein), little is known about its analytic properties, the major 
work having been done in developing algorithms for numerical calculations [Fet74]. 
Dan Lozier, remarked the need for new techniques in the computation of inverf (z) 
[Loz96]. 

In this paper, we analyze the asymptotic behavior of the derivatives inverf (z) 
for large values of n, using a discrete WKB method [CC96]. In Section 2 we present 
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some properties of the derivatives of inverf (z) and review our previous work on 
nested derivatives. In Section 3 we study a family of polynomials Pn{x) associated 
with the derivatives of inverf {z), which were introduced by L. Carlitz in [Car63]. 
Theorem 3.3 contains our main result on the asymptotic analysis of Pn{x)- In Sec- 
tion 4 we give asymptotic approximations for inverf (-j) and some numerical 
results testing the accuracy of our formulas. 

2. Derivatives 

Let us denote the function inverf (z) by 3{z) and its derivatives by 



(2.1) 



dn = — inverf (z) 

dz'^ ^ ^ 



n ■■ 



0,1,.... 



z=0 



Since erf(z) tends to ±1 as ^; — > ±oo, it is clear that inverf (z) is defined in the 
interval (—1, 1) and has singularities at the end points. 

Proposition 2.1. The function 3{z) satisfies the nonlinear differential equa- 
tion 

(2.2) 3" -23{3'f ^0 

with initial conditions 

3(0) =0, 3'{0) = ^. 
Proof. It is clear that 3(0) = 0, since erf(O) = 0. Using the chain rule, we 



(2.3) 

F 
have 



_L_.:^expP^[erfWl} 



and therefore 

(2.4) 



a'=^exp(3^). 



Setting z = we get 3'{0) = and taking the logarithmic derivative of (2.4) the 
result follows. □ 

To compute higher derivatives of 3(z), we begin by establishing the following 
corollary. 

Corollary 2.2. The function 3{z) satisfies the nonlinear differential-integral 

equation 



(2.5) 



Z 

3'iz) J 3{t)dt = -l + ^3'{z). 



Proof. Rewriting (2.2) as 



and integrating, we get 



3 = 



1 3" 



2(30 



Z 

/ 



3{t)dt = - 



1 



1 



3'{z) 3'{0) 



1 2 
+ 



3'{z) V^, 



and multiplying by 3'{z) we obtain (2.5). 



□ 
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Proposition 2.3. The derivatives of3{z) satisfy the nonlinear recurrence 

n-l 



(2.6) 



dn 



+1 



k=0 



n 
k + 1 



dkdn-k, n = 1,2, . . . 



with do = and di = 
Proof. Using 

and di = in (2.5), we have 



z 



n=l 



°° Z"' 1 
/ , dn-l—r 1= 



.n=l 



or 



n=l n=2 



E(, , 1 rffcd„_fc -j- ^Vdn+i— =0. 
V A; + 1/ n! vtt n! 

Comparing powers of 2;", we get 



, — It — z, 



or 



E 

fe=0 



fe=0 

n 
A;+ 1 



A; + 1 



dkdn-k l=dn+l — 



dkdn-k 7=C^n+l = 0. 



□ 



Although one could use (2.6) to compute the higher derivatives of inverf (z) , 
the nonlinearity of the recurrence makes it hard to analyze the asymptotic behavior 
of d„ as ri ^ cx). Instead, wc shall use an alternative technique that we developed 
in [DomOSa] and we called the method of "nested derivatives". The following 
theorem contains the main result presented in [DomOSa] . 



Theorem 2.4. Let 

H{x)=h-\x), f{x) 

Then, 



h'{x) 



zo = h{xo), [/(xo)! G (0,00) . 



Hiz)=xo + f{xo)J2^''~'if] M 



{z - zoT 



where we define (x), the n^^ nested derivative of the function f{x), by (x) 



(^) = ^ im X i^)] , n = 0, 1, . . . . 



1 

(2.7) 



The following proposition makes the computation of D" [/] (xo) easier in some 
cases. 
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Proposition 2.5. Let 

(2.8) s)"[/](^) = f;A^l^_pl!, /(^) = gB,i^^p^. 

fe=0 ■ k=0 

Then, 

k+l 

(2.9) A-,+^ = {k + l)J2M+i-jBj. 
Proof. Prom (2.8) we have 

(2.10) /(^)s"[/](^)=f:a^^^^p^, 

fe=0 

with 

k 

(2.11) a^ = ^A^_,i?,. 
Using (2.8) and (2.10) in (2.7), we obtain 

oo , oo oo 

/s=0 /s=0 fe=0 

and the result follows from (2.11). □ 

To obtain a linear relation between successive nested derivatives, we start by 
establishing the following lemma. 

Lemma 2.6. Let 
Then, 

(2.13) gn+i{x)=g'^ix) + {n + l)Q^gn{x), n = 0,l,.... 

Proof. Using (2.7) in (2.12), we have 

_ (^) _ ^[/(a^)xg»[/](x)] 



[i?n {x) r+' {x)] _ g'^ {x) (x) + gn (x) (n + 1)/" (x) fix) 



fn+1 (a.) (^) 

and the result follows. □ 
Corollary 2.7. Let 

H{x) = h-\x), /(xO-^, zo = /i(xo), |/(a;o)| e (0,oo). 

Then, 

(2.14) ^(^o) = [/(xo)r5n-i(^o), n=l,2,.... 
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For the function h{x) = crf(z), wc have 

(2.15) /(^) = /7M = '^^^^^^^^' 

and setting a;o = we obtain zo = erf(O) = 0. Using the Taylor series 

in (2.9), we get 



3=0 •' 

with Al defined in (2.8). Using (2.15) in (2.13), we have 

(2.16) gn+i{x)=g'n{x)+2(n+l)xgn{x), n = 0, 1,..., 
while (2.14) gives 

(2.17) dn=(^\ Qn-m, n = l,2,.... 



In the next section wc shall find an asymptotic approximation for a family of 
polynomials closely related to Qn {x). 

3. The polynomials P„{x) 

We define the polynomials Pn{x) by Pq{x) = 1 and 

(3.1) P„(x)=5n(^)2-t. 

(3.2) Pn+i{x)^P^{x) + {n+l)xPn{x), 
The first few P„ [x) are 

Pi{x) ^ X, P2{x) = 1 + 2a;^ P^ix) = 7,x + 6a;^ .... 
The following propositions describe some properties of P„ (x) . 
Proposition 3.1. Let 



(3.3) P„(x) = E^fe^"" 



-2k 
k=0 

where [■} denotes the integer part function. Then, 

(3.4) Co" = n! 
and 

n-l jfc-1 i2-l k ■ _9- I 9 

(3.5) C," = n!E E •■•En^^-ff^' ^=1-- 

jfc=Ojib-i=0 ji=Oi=l .^'"^ 
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Proof. Using (3.3) in (3.2) we have 

^„+i^„+i-2fc^ (n - 2fc) + Yl + 1) Cfc 

0<2fe<n+l 0<2k<n 0<2fc<n 

= Yl C-^-i - 2fc + 2) + Y (n+l)Cfc2;"+^"^*^- 

2<2fc<n+2 0<2fe<n 

Comparing coefficients in the equation above, we get 
(3.6) 



2fc 



Cn+1 /^n 
— ' 



Cr' = («-2fc + 2)C^_i + (n+l)C,", k = l,..., 



(3.7) ^"+^ 
and for n = 2m — 1, 

^2m _ ^2m-l _ i 9 

'-'m — '-'m-1 ' TO— i,/, 

Prom (3.6) we immediately conclude that Cg = n!, while (3.7) gives 

n-l 

(3.8) C^=nlY 



i - 2fc + 2 • 



Setting A; = 1 in (3.8) and using (3.4), we have 

n— 1 . n— 1 

(3-9) ^" = -'E7-TTj!^o=-!E 7TT- 



Similarly, setting A: = 2 in (3.8) and using (3.9), we get 



n-i . ^ r j-i 



J - 2 I 



i=o ^-^ 'I i=0 J j=0 i=0 



and continuing this way wc obtain (3.5). 



□ 



Proposition 3.2. The zeros of the polynomials Pn{x) are purely imaginary 
for n > 1. 

Proof. For n = 1 the result is obviously true. Assuming that it is true for n 
and that Pn{x) is written in the form 



(3.10) 



Pn{x) =nlY[iz- Zk), Re{zk)=0, l<k<n, 



we have two possibilities for z*, with P„+i(^*) = 0: 

(1) z* = Zk, for some 1 < A: < n. 

In this case, Re(^!*) = and the proposition is proved. 

(2) z* ^ zu, for alll < fc < n. 

From (3.2) and (3.10) we get 

= i m [Pnix)] + (n + 1). ^ t + + 



k=l 



Pn{x) dx 

Evaluating at 2; = 2*, we obtain 

n 

= Y- + {n + l)z* 



k=l 
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and taking Re(»), we have 
= Re 



f +{n+l)Re{z*) = Re{z*) 

\Z* - Zk\ 



E 

k=l 



which impUes that Re{z*) = 0. 



" 1 

I ^* 



+ n+ 1 



- Zk\ 



□ 



3.1. Asymptotic analysis of Pn{x). We first consider solutions to (3.2) of 
the form 

(3.11) P„(a;)=nU("+i)(a;), 

with a; > 0. Replacing (3.11) in (3.2) and simplifying the resulting expression, we 
obtain 

A^ix) =A'{x)+xA{x), 

with solution 



(3.12) 



A{x) = exp ( -y 



C 



-, -1 



for some constant C. Note that (3.11) is not an exact solution of (3.2), since it docs 
not satisfy the initial condition Po{x) = 1. To determine C in(3.12), we observe 
from (3.4) that 

(3.13) Pn{x) ~ n!a;", x ^ oo. 

As a; ^ oo, we get from (3.12) 



ln[A(a;)]~---ln(C- 



exp 



a; ^ oo, 



(c-yf)x' 

which is inconsistent with (3.13) unless C = y^. In this case, we have 



(3.14) 

matching (3.13). Thus, 

(3.15) A{x) = 

Since (3.11) and (3.14) give 



Aix) ~ X H — , a; — > oo, 

X 



2 / x2 
exp 



1 — erf 



Pn{x) ^ n\x 



a; ^ 00, 



instead of (3.13), we need to consider 

(3.16) Pn{x) = n\A^'^+^\x)B{x,n). 

Replacing (3.16) in (3.2) and simplifying, we get 

1 dB 



B{x,n + l) = B{x,n) + 



A{x){n + l) dx 



{x,n). 
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Using the approximation 



0J3 1 &^ B 

B{x,n+l) = B{x,n) + + 2"^(^'"') + ' 



we obtain 



1 



dB 



whose solution is 
(3.17) 



as 

'dn ~ A{x){n + l)'dx' 
B{x,n)=F ""^^ 



for some function F{u). Matching (3.16) with (3.13) requires 

1 



(3.18) 



B{x, n) 



a; — > oo. 



Since in the limit as a; — > oo, with n fixed we have 



In 



n + 1 



X 

y 



(3.17)-(3.18) imply 

Therefore, for a; > 0, 

(3.19) 

with 

$ {x, n) = 



F{u) = . 
Pn{x) {x,n) , n — > oo, 



2 cxp 
1 - erf 



(*) 



n + l 



2 In 



n+l 



1 - erf 



(^) 



Prom (3.3) we know that the polynomials Pn{x) satisfy the reflection formula 

(3.20) P„(-a;) = (-l)"P„(a;). 

Using (3.20), we can extend (3.19) to the whole real line and write 

(3.21) P„(a;) ~n![$(a;,n) + (-l)"$(-a;,n)], n ^ oo. 

In Figure 1 we compare the values of Pio(a;) with the asymptotic approximation 

(3.21) . 

We see that the approximation is very good, even for small values of n. We 
summarize our results of this section in the following theorem. 

Theorem 3.3. Let the polynomials Pn{x) be defined by 

Pn+l{x) = P^{X) + (n + 1) xPn{x), 

with Pq{x) = 1. Then, we have 

(3.22) P„(a;)~P„(a;)~n![$(a;,n) + (-l)"$(-a;,n)], n ^ oo. 
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where 



(3.23) 



$ {x,n) = 



2 exp 

TT 



n+1 



2 In 



4. Higher derivatives of inverf (z) 
From (2.17) and (3.1), it follows that 



(4.1) 



P„_i(0), n=l,2,..., 



where (i„ was defined in (2.1). Using Theorem 3.3 in (4.1), we have 



dn 



1 



\/2 VV 2 
as n — > oo. Using (3.23), we obtain 

ri! 2n-^ln(n 



$(0,n-l) l + (-l) 



(4.2) 



i+(-ir 



n — > oo. 
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Setting n = 2N + 1 in (4.2), we have 

{2N + 1)! ^ (2iV + 1) VR2iV + 1) ' 



4.1. Numerical results. In this section we demonstrate the accuracy of the 
approximation (4.2) and construct a high order Taylor scries for inverf (x) . In Figure 

2 we compare the logarithm of the exact values of inverf (x) and our 

x—O 

asymptotic formula (4.2). We see that there is a very good agreement, even for 
moderate values of n. 

Using (2.6), we compute the exact values 

1 1 , 1 3 , 7 5 , 127 7 , 4369 9 
and form the polynomial Taylor approximation 

4 2fc+l 
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o / 

/ 

° 1 


\ *° 


o- / 

d/ 

i 0.8- 


\b 
e 


0.6 




0.4 




0.2 





-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

X 



Figure 3. A sketch of 3M.^ (solid curve), (+++) 
and (ooo). 

inverf(a;) ' 



In Figure 3 we graph 3^.,- and ZiMl^ivM f ^ ^ ;^o, 20, where 

o o iTiverT I .T ) invert ( x i ' ' ' 



inverf (x) 
N 



(4.4) 



i?^(a;) = ^ 



inverf(x) 

,2/c+l 



^5(2iV + l) Vln(2iV + l) 



A/' = 5,6,. 



The functions are virtually identical in most of the interval (—1,1) except for 
values close to x = ±1. We show the differences in detail in Figure 4. Clearly, the 
additional terms in R2o{x) give a far better approximation for a; ~ 1. 

In the table below we compute the exact value of and optimal asymptotic 
approximation to inverf (x) for some x: 



X 


inverf (.x) 


mx)+RN{x) 


N 


0.7 


.732869 


.732751 


6 


0.8 


.906194 


.905545 


7 


0.9 


1.16309 


1.16274 


11 


0.99 


1.82139 


1.82121 


57 


0.999 


2.32675 


2.32676 


423 


0.9999 


2.751UG 


2.73103 


3G85 
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0.4- 



0.2 



0.9 0.92 0.94 0.96 0.98 1 

X 

Figure 4. A sketch of 3M.^ (solid curve), (+++) 
and (ooo). 

inverf(a;) ' 



Clearly, (4.4) is still valid for a; — > 1, but at the cost of having to compute many 
terms in the sum. In this region it is better to use the formula [DomOSb] 



inverf {x) 



\ 



■LW 



TT (a; — 1) 



where LW(-) denotes the Lambert-W function [CGH+96], which satisfies 

LW(a;) cxp [LW(a;)] = x. 
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